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ABSTRACT 

We present simultaneous observations by EUVE, ASCA, and RXTE of the type 1 Seyfert galaxy NGC 5548. 
These data indicate that variations in the EUV emission (at ~ 0.2keV) appear to lead similar modulations in 
higher energy (;> 1 keV) X-rays by ^10-30ks. This is contrary to popular models which attribute the correlated 
variability of the EUV, UV and optical emission in type 1 Seyferts to reprocessing of higher energy radiation. 
This behavior instead suggests that the variability of the optical through EUV emission is an important driver for 
the variability of the harder X-rays which are likely produced by thermal Comptonization. We also investigate 
the spectral characteristics of the fluorescent iron Ka line and Compton reflection emission. In contrast to prior 
measurements of these spectral features, we find that the iron Ka line has a relatively small equivalent width 
{Wko ~ 100eV) and that the reflection component is consistent with a covering factor which is significantly less 
than unity (H/2tt ~ 0.4-0.5). Notably, although the 2-10 keV X-ray flux varies by ~ ±25% and the derived 
reflection fraction appears to be constant throughout our observations, the flux in the Fe Ka line is also constant. 
This behavior is difficult to reconcile in the context of standard Compton reflection models. 

Subject headings: galaxies: individual (NGC 5548) — galaxies: Seyfert — X-rays: galaxies 



1. INTRODUCTION 

Over the past decade, a moderately coherent picture of 
the broad band continuum emission from Seyfert AGNs has 
emerged. In this model, an accretion disk around a super- 
massive black hole (M ~ 10 7 -10 8 M Q ) produces thermal emis- 
sion primarily in the optical and ultraviolet; and a hot, Comp- 
tonizing corona above the disk up-scatters these photons to pro- 
duce X-rays with energies ~ 1-100 keV. Furthermore, the exis- 
tence of strong fluorescent iron Ka lines at ~ 6.4keV and the 
so-called "Compton reflection humps" above ~ lOkeV in the 
spectra of type 1 Seyferts indicate that the disk, or some other 
cold, optically thick material reprocesses the hard X-rays. 

Central uncertainties in this model are the geometry of the 
corona and the disk and how the radiation from one compo- 
nent affects the properties and emission of the other. A cold 
(T <J 10 5 K), optically thick disk will absorb a large fraction 
of the incident X-rays and will likely be significantly heated, 
thus reprocessing the absorbed energy as additional thermal 
photons. These photons will be available for Compton up- 
scattering and may provide additional cooling of the coronal 
electrons thereby decreasing the temperature of the Comptoniz- 
ing medium. Hence, the size, shape, and location of the corona 
relative to the disk will affect the temperature and spectral prop- 
erties of both components. 

A related difficulty is the origin of the soft X-ray excess at en- 
ergies ;> 0. 1 keV which has no plausible explanation within the 
context of standard thin disk models (Koratkar & Blaes 1999). 
Any reasonable range of disk temperatures for sub-Eddington 
accretion is too low to account for this emission as an exten- 
sion of the disk thermal emission, and an extrapolation of the 



standard hard X-ray thermal Comptonization spectrum down to 
these energies severely underpredicts the observed flux. How- 
ever, Magdziarz et al. (1998) have been able to model the soft 
X-ray excess and most of the optical/UV as emission from very 
Thomson thick (r ~ 30), warm (200 eV) Comptonizing clouds 
above the disk. Despite the ad hoc nature of this model, it 
makes the definite prediction that the optical, UV and EUV flux 
should be highly correlated. 

Correlated variability studies between different wavebands 
have been the traditional means of determining which spectral 
components provide the impetus for others. NGC 5548, in par- 
ticular, has been involved in several multiwavelength monitor- 
ing campaigns, the best known of which were efforts to use 
optical/UV spectra to characterize the Broad Line Region via 
reverberation mapping (Clavel et al. 1991; Korista et al. 1995). 
As a by-product of these investigations, Peterson et al. (1991) 
compared the light curves of contemporaneous optical (4870A) 
and UV (1350A) data from the earlier of these campaigns and 
showed that essentially no lag exists between spectral bands 
within the temporal resolution of the observations, ±2 days (see 
also Clavel et al. 1991). If the observed optical/UV flux in 
type 1 Seyferts is thermal emission due to the local conversion 
of gravitational energy in thin accretion disks, then the charac- 
teristic radii from which continuum emission at these two wave- 
lengths are produced will be quite different. Any correlated 
variability between these two bands would then be limited by 
the sound crossing time in the disk flow. However, Courvoisier 
& Clavel (1991) have shown that the upper limits on the opti- 
cal/UV lag for NGC 5548, as well as for Seyferts NGC 4151 
and Fairall 9, are several orders of magnitude smaller than the 
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sound crossing times. 

As an explanation for this near simultaneity of the opti- 
cal and UV, Clavel et al. (1992; see also Malkan 1991) sug- 
gested that the variable emission is due to reprocessing of X- 
rays from a source located very near the central black hole. 
This interpretation would also be consistent with the obser- 
vations of NGC 5548 by the EUVE satellite during the 1993 
/[/£/ffi>77ground-based BLR reverberation mapping campaign 
(Marshall et al. 1997; Korista et al. 1995). Marshall et al. ob- 
tained contemporaneous EUVE data and found that the EUV 
(at ~ 0.16keV) also varies with the optical and UV to within 
<J 0.25 days. This result links the soft X-ray excess with the 
lower wavebands and supports either a reprocessing model or 
one such as that of Magdziarz et al. (1998). In conflict with 
the former scenario are the contemporaneous observations of 
the type 1 Seyfert NGC 7469 by RXTE and IUE (Nandra et al. 
1998). Over a 30 day baseline, Nandra et al. found that no sim- 
ple relationship exists between the variability in the X-rays and 
UV. Light curves in both wavebands exhibit large, factor of ~ 2 
modulations on 10 day time scales, but the peaks in the UV ap- 
pear to lead the corresponding maxima in the X-rays by ~ 4 
days while the troughs in the two wavebands are nearly simul- 
taneous. Nandra et al. suggest that this behavior may be due 
to several reprocessing regions which interact at different time 
scales depending on the flux state. Previous correlated UV/hard 
X-ray observations of NGC 5548 by Ginga and IUE were only 
able to place an upper limit of ~ 6 days for the characteristic lag 
between these two energy bands (Clavel et al. 1992). An ear- 
lier observation by EXOSAT did show an apparent lag of ~ 5 ks 
between the soft (0.05-2.5 keV) and hard (2-10 keV) X-ray 
bands, but this was only seen during one of 12 observations, 
albeit the longest one of 60 ks, and with fairly poor statistics 
(Kaastra & Barr 1989; see also Nandra et al. 1991). 

Another aspect of the behavior of NGC 5548 which bears 
upon the disk/corona geometry is the putative correlation be- 
tween the derived reflection fraction and the spectral index 
of the underlying hard continuum (Magdziarz et al. 1998; 
Zdziarski, Lubinski, & Smith 1999). Spectral analyses of X- 
ray data from Galactic black hole candidates and Seyfert AGNs, 
including NGC 5548, show that the fitted reflection fraction is 
systematically larger for softer underlying spectra. The natural 
interpretation of such a correlation is that the reflecting medium 
is an important source of soft seed photons which affect the 
cooling of the Comptonizing medium. A larger covering factor 
implies more soft photons, causing a lower temperature coronal 
region and therefore a softer spectrum. Since the fluorescent 
iron line emission presumably also originates in this same cold, 
Compton reflecting material, it provides an additional diagnos- 
tic to explore the relationship between the covering factor and 
the spectral hardness. 

In order to investigate the roles of the various components 
in the Seyfert disk/corona system of NGC 5548 further, we 
have conducted simultaneous EUVE, ASCA, and RXTE mon- 
itoring observations. This combination of instruments allowed 
us to study the soft X-ray excess, its relationship to the harder, 
Comptonized X-ray emission, and the higher energy signatures 
of X-ray reprocessing in accretion disks, the fluorescent iron 
Ka line and the Compton reflection component. We describe 
the observations by the various instruments in § 2. In § 3, 
we present an analysis of the correlated variability analysis be- 
tween different energy bands, and in § 4 we examine the power 
spectra of the various bands. In § 5, we analyze the X-ray spec- 
tral properties in detail during each observation epoch. We dis- 



cuss the implications of our observations for disk/corona mod- 
els in § 6 and present our summary and conclusions in § 7. 

2. THE OBSERVATIONS AND DATA REDUCTION 

Great effort was made to ensure that the observations by the 
various telescopes in this campaign were as simultaneous as 
possible. Unfortunately, because of differences in time allo- 
cation, conflicts with other scheduled observations, differing 
visibility constraints and other unforeseen complications, ab- 
solute simultaneity over our proposed 23 day baseline was not 
possible. In particular, several data gaps of a few x lOks exist 
in the RXTE data because of scheduling conflicts; the EUVE 
spacecraft went into "safepoint" mode just prior to one of our 
shorter monitoring observations; and the final RXTE observa- 
tion which took place over a 200 ks time period had no con- 
temporaneous EUVE or ASCA monitoring because of visibility 
constraints. Nonetheless, these data represent some of the high- 
est time resolution broad band monitoring of any AGN. 

In order to provide uniformity in referring to the various 
pointings performed by the three satellites, we have adopted 
a numbering scheme based on the originally scheduled simul- 
taneous observations. In this scheme, the June 15, June 18- 
24, July 1 and July 7 observations are numbered 1-4, with the 
June 18-23 (obs. 2) subdivided as we discuss below. We refer 
to all the other observations taken by individual instruments by 
the dates on which they took place. 

2.1. EUVE 

Our strategy for using EUVE to monitor the soft X-ray flux 
from NGC 5548 was identical to that employed by Marshall et 
al. (1997). Photometry was carried out using the Deep Survey 
(DS) detector with the Lex/B filter, and the pointing was off- 
set from nominal by 0.3 degrees along the direction of the SW 
spectrometer dispersion axis. This was done in order to avoid 
the deadspot in the DS detector while largely preserving the ca- 
pabilities of the SW spectrometer. In extracting the count rates, 
we employed the analysis of Marshall et al. which accounts for 
dead-time and "Primbsch" corrections as in the standard EUVE 
analysis procedure and which also applies a vignetting correc- 
tion necessitated by the offset pointing. Another important fea- 
ture of this analysis is its implementation of a maximum likeli- 
hood technique using the instrumental point spread function to 
determine the source and background counts rather than speci- 
fying source and background apertures. Further information on 
these procedures can be found in Marshall et al. (1997). 

Table [j] gives the details of our EUVE-DS observations. The 
spectral densities we report account for a Galactic column den- 
sity of Nh — 1 .7 x 10 20 cm" 2 and assume a power-law spectrum 
with index 0-2 (/„ oc v~ a ). For the narrow effective bandpass 
we consider (70-90A), the resulting spectral densities we com- 
pute are not very dependent on the value of a. The 1998 June 
2, June 9, and July 1 observations were simultaneous with Lick 
observations scheduled on those dates, although the June 2 Lick 
observation did not occur due to bad weather. As we noted ear- 
lier, the EUVE spacecraft went into "safepoint" mode just prior 
to our June 15 joint ASCA/ RXTE observations (cf. Tables || 
& [|), and we did not obtain simultaneous EUVE data for our 
July 7 ASCAIRXTE observations because of a scheduling error. 

2.2. ASCA 
2.2.1. Basic data reduction 
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ASCA observed NGC 5548 on five occasions during this 
campaign. The dates, good exposure times, count rates and 2- 
lOkeV fluxes are reported in Table ||. ASCA possesses four de- 
tectors; two solid-state imaging spectrometers (SISO and SIS 1) 
and two gas-imaging spectrometers (GIS2 and GIS3). These 
detectors are located at the focal points of four independent but 
coaligned X-ray telescopes. During our campaign, data from 
all four detectors were obtained. SIS data taken in both BRIGHT 
and FAINT mode were combined in order to maximize the to- 
tal signal and a standard GRADE selection was performed in 
order to reduce the effects of particle and instrumental back- 
ground. Data from the SIS were further cleaned in order to 
remove the effects of hot and flickering pixels and subjected to 
the following data-selection criteria: the satellite should not be 
in the South Atlantic Anomaly (SAA), the object should be at 
least 5° above the Earth's limb, the object should be at least 25° 
above the day-time Earth limb, and the local geomagnetic cut- 
off rigidity (COR) should be greater than 6 Ge V/c. Data from 
the GIS were cleaned to remove the particle background and 
subjected to the following data-selection criteria: the satellite 
should not be in the SAA, the object should be at least 7° above 
the Earth's limb and the COR should be greater than 7 GeV/c. 
These steps were performed using XSELECT vl .4 which is part 
of the ftools v4. 1 package. 

Images, light curves and spectra were extracted from the 
good data using a circular region centered on the source. For 
the SIS, an extraction radius of 3 arcmins is used whereas a ra- 
dius of 4 arcmins is used for the GIS. These regions are suffi- 
ciently large to contain all but a negligible portion of the source 
counts. Background spectra were extracted from source free 
regions of the same field of view within each of the four ASCA 
instruments. Background regions for the SIS were taken to be 
rectangular regions along the edges of the source chip whereas 
annular regions were used to extract GIS background. Approx- 
imately 1 % of the photons in the source extraction region were 
found to be from the background. 

2.2.2. The effects of RDD 

The effects of radiation damage to the CCD's of both SIS 
detectors have reduced the ability of the on-board software to 
estimate and correct for dark current. This problem is known 
as residual dark distribution (RDD) and, in these late stages of 
the ASCA mission, is becoming a major concern. Data taken in 
FAINT mode can have the effects of RDD partially corrected for 
during the ground-based reduction process. Unfortunately, due 
to telemetry constraints, the majority of our data were taken in 
BRIGHT mode and cannot be corrected. 

Operationally, RDD renders the soft end of the SIS spectrum 
(< 1 keV) untrustworthy. Comparing our SISO and SIS1 spec- 
tra reveals major discrepancies between these two instruments 
below 1 keV which are readily attributed to RDD effects. For 
this reason, we restrict ourselves to the 2-10 keV band when 
performing spectral fitting with SIS data. 

2.2.3. An SIS/GIS discrepancy 

The spectral fitting described below was initially performed 
using the 2-10keV data from all four ASCA instruments. How- 
ever, the GIS spectra were found to strongly disagree with the 
SIS spectra at high-energies (8-10keV) — the two GIS spectra 
consistently lay above the two SIS spectra in each and every ob- 
servation of this campaign. When these spectra are fitted with a 
spectral model consisting of an absorbed power-law model and 
a Gaussian emission line (see below), residuals are present in 



the GIS spectra at energies coincident with the Gold M-edge 
(2.2 keV) and Xenon L-edge (4.8 keV). This strongly indicates 
the use of an incorrect GIS gain factor. No linear gain correc- 
tion could eliminate these residuals simultaneously. Faced with 
the conclusion that there is a non-linear gain problem in our 
GIS data, we choose not to use GIS data in any of the spectral 
fitting reported below. 

Further support for the hypothesis that the GIS is in error 
is given when RXTE-PCA data are fitted simultaneously with 
these ASCA data. In the overlap band (4-10keV), the PCA 
spectra generally agree well with SIS spectra and disagree with 
GIS spectra. However, to our knowledge, the ASCA and RXTE 
instruments have not been cross-calibrated. As we discuss be- 
low, we have therefore allowed the relative normalization for 
the spectral models which are fit to data from each instrument 
to vary, and we have also allowed the spectral index of the hard 
continuum to vary separately for the SIS and PCA data. 

2.3. RXTE 

Standard extraction methods were employed to determine 
the PCA and HEXTE light curves and spectra. The PCA in- 
strument consists of five collimated (1° FWHM) proportional 
counter units (PCUs) which are numbered 0^1. Each PCU con- 
tains three multianode detector layers with a mixture of xenon 
and methane gas, has a bandpass of 2-60 keV and a geomet- 
ric collecting area of ~ 1400cm 2 . The energy resolution of 
each detector is ~ 8% FWHM at 6.6 keV (Glasser, Odell, & 
Seufert 1994). PCA data were collected only from PCUs 0- 
2, since PCUs 3 and 4 were turned on for a smaller fraction 
of the on-source time due to breakdown. This source is rel- 
atively weak, so in order to maximize signal-to-noise, we ac- 
cumulated events only from the top xenon/methane layer of 
the PCA. As recommended by the PCA team, data were dis- 
carded for the 30 minutes following a SAA passage, during 
Earth occultation (i.e., when the Earth elevation angle is < 10°), 
and when there is severe electron contamination. The vari- 
ability and spectral analyses were greatly facilitated by the re- 
cently developed faint source (< 40cps) background model. It 
is implemented in the PCABACKEST program and is contained 
in the model files pca_bkgd_f aintl7_e03v03 .mdl and 
pca_bkgd_faint240_e03v03 . mdl. The first file con- 
tains information on the variation on the "L7" rate which is a 
background rate made up of a combination of 7 rates in pairs 
of adjacent signal anodes, while the second file is related to the 
integrated recent particle doses during SAA passages, as mea- 
sured by the HEXTE particle monitor. SAA passages and doses 
are recorded in pca_saa_history_v8. Response matrices 
were constructed using PCARMF v3.5 and PCARSP v2.36. For 
spectral fitting, we restricted the energy range of the PCA data 
from 3 to 20 keV 

The HEXTE instrument consists of two clusters of four 
Nal/CsI-phoswich scintillation counters which are sensitive 
from 15 to 250 keV with an open area of ~ 900 cm 2 for each 
cluster. Beam switching or "rocking" of these clusters be- 
tween source and background fields provided direct measure- 
ments of the HEXTE background with an on-source duty cycle 
of ~ 60%. As recommended by the HEXTE team, data were 
discarded for the 2 minutes following a SAA passage and dur- 
ing Earth occultation. Standard deadtime corrections were ap- 
plied to the Science Archive data. In order to maximize signal- 
to-noise in the light curves, only counts from absolute detec- 
tor channels 30-122 were used. As can be seen from Table |[ 
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the HEXTE count rates were very low, < 1 cps throughout our 
observations. Hence, statistically meaningful HEXTE spectra 
could only be accumulated over relatively long integrations, 
2> 4 x 10 5 s, which is substantially longer than the variability 
time scale seen in the PCA data. 

2.3.1. IE 1415.6+2557 

A minor complication in the RXTE analysis is the existence 
of a contaminating source within the field-of-view of the PCA 
and HEXTE collimators, 0.5° away from NGC 5548. The 
source is a BL Lac object, IE 1415.6+2557, and has been pre- 
viously observed by ROSAT (Nandra et al. 1993). As an addi- 
tional precaution, since Nandra et al. reported several additional 
sources in the ROSAT-PSPC FOV surrounding NGC 5548, we 
performed a PCA scanning observation of the field to deter- 
mine if any other contaminants were present (cf. Marshall et al. 
1998). Using the collimator response model of the PCA which 
was provided to us by the PCA instrument team, we fit the scan 
observations and found that we could account for all the counts 
in the field adequately with just two sources: NGC 5548 and 
IE 1415.6+2557. 

During the first four of our RXTE observing periods, the con- 
tribution by IE 1415.6+2557 to our NGC 5548 measurements 
was determined from ~ 3-5 ks pointings directed at IE 1415.6 
+2557 just prior to and after each NGC 5548 observation pe- 
riod. This bracketing was done in case IE 1415.6+2557 exhib- 
ited any variability on 30 ks time scales. We used the afore- 
mentioned collimator model to fit the count rates from these 
two sources simultaneously during the beginning and ending 
sub-portions of each observation. We then extracted spec- 
tra from these sub-portions for each source. The spectra for 
IE 1415.6+2557 were fit with XSPEC vlO.00 using the associ- 
ated background-subtracted NGC 5548 spectrum as a "correc- 
tion file", appropriately scaled to yield the fitted count rates. 
Table |] gives the estimated count rates and spectral fits for 
IE 1415.6+2557 obtained from this procedure. When there are 
significant differences in the count rates of IE 1415.6+2557 for 
the two bracketing pointings, they are reported separately, i.e., 
as 2a/b and 3a/b. 

For the observations beginning 1998 August 16, we al- 
tered our strategy for estimating the flux and spectrum from 
IE 1415.6+2557. Rather than relying on scaled spectral correc- 
tion files extracted from direct pointings at NGC 5548 which 
are in turn contaminated by flux from IE 1415.6+2557, we took 
"background" measurements at the symmetric location on the 
opposite side of NGC 5548. For the same roll angles, these 
pointings should contain the same contribution from NGC 5548 
as the direct IE 1415.6+2557 pointings, but since they are ~ 1° 
away, they should be uncontaminated by IE 1415.6+2557 it- 
self. Correction files from these pointings were then used in the 
spectral analysis. We find that the spectrum of IE 1415.6+2557 
during the August 16 observations was somewhat harder and 
the count rate in the PCA had increased by ~ 40%. We note 
that this increase is not due to our altered observing strategy as 
we have also applied the collimator model fitting method and 
confirm this result. The collimator response for pointings di- 
rected at NGC 5548 to photons from IE 1415.6+2557 is ~ 0.33, 
so the contaminating contribution to the NGC 5548 light curves 
will be ~ 2 cps with a modulation of <J 0.8 cps, while the PCA 
count rates for NGC 5548 are typically 20-25 cps. Therefore, 
in our variability analysis of NGC 5548, we can safely ignore 
contributions from IE 1415.6+2557. 
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FIG. 1 .— EUVE, ASCA-SIS, and RXTE-PCA light curves for the long 
observation. All light curves have been normalized to unity mean. For 
clarity, the ASCA and RXTE curves were then offset by 0.75 and 1.5 re- 
spectively. The dashed lines show the boundaries of the four intervals 
which were used for studying the intra-observation spectral variability. 

3. CORRELATED VARIABILITY 

In Fig. H, we show the EUVE, ASCA, and RXTE-PCA count 
rate light curves for the longest of our contemporaneous obser- 
vation periods. The most distinct feature in all three of these 
light curves is the large step at ~ 2 x 10 4 s in the figure. This 
step is more pronounced (~ 40%) at EUV energies than in the 
harder X-rays (~ 20% for the RXTE-PCA data). After this step, 
in theRXTE light curve there are several local maxima which re- 
cur on ~ 50 ks time scales; similar modulations also appear in 
the EUVE light curve, but with less certainty due to the rela- 
tively poor statistics. As a first order analysis of the spectral 
variability across these wavebands, we have computed cross- 
correlation functions (CCFs) for each pair-wise combination. 

The method we use is the Z-transformed Discrete Cross- 
Correlation Function (ZDCF) of Alexander (1997) based upon 
the DCF method of Edelson & Krolik (1988). In the latter pro- 
cedure, which was designed to accommodate unevenly and dif- 
ferently sampled data trains, pair-wise combinations of mea- 
sured flux values from the two light curves are binned accord- 
ing to their relative delays or lags. For the data pairs within 
each lag bin, a quantity which is essentially Pearson's linear 
correlation coefficient is computed (Press et al. 1992). In the 
continuum limit, this procedure is equivalent to the standard 
definition of the CCF Alexander (1997) proposed two modifi- 
cations of this method: first, the lag bins are defined to contain 
a specific number of data pairs which should be > 1 1 in order to 
ensure convergence; second, the errors are estimated by apply- 
ing the Z-transform to the correlation function values and using 
the number of pairs for the given bin and known properties of 
the Z-transform to estimate the errors. For identical binning, the 
ZDCF estimates of the CCF are equal to those produced by the 
DCF, but the uncertainties are better behaved. Auto-correlation 
functions can also be computed by this method. 

Figure | shows the ZDCFs for EUVE vs 0.5-1 keV ASCA- 
SIS, EUVE vs RXTE-PCA, and 0.5-1 keV ASCA-SIS vs RXTE- 
PCA. In each case, a positive time delay indicates the latter light 
curve lagging the former. In order to obtain a quantitative esti- 
mate of the relative lags, we fit the ZDCF values in the vicinity 
of the peak with a parabola (solid curves in Fig. ^). This func- 
tional form has no specific physical meaning — it is merely a 
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FIG. 2. — Cross-correlation functions computed using the ZDCF 
(Alexander 1997). The bands which are compared are EUVE (0.14- 
0.18 keV), ASCA (0.5-1 keV), and RXTE-PCA (2-20 keV) (cf. Fig. |l|). 
The solid curves are the parabolas which were fit to find the location of 
the characteristic lags. 



convenient parametrization of the ZDCF shape near the max- 
imum. From these fits, we see that higher energy light curve 
modulations are delayed relative to lower energy ones with the 
variations in the 2-20 keV RXTE-PCA data lagging those of 
the EUVE by - 35 ks and those of the 0.5-1 keV ASCA data 
by ~ 5ks. Similarly, the 0.5-1 keV ASCA data lag the EUVE 
by ~ 13 ks. We determine confidence limits using Monte Carlo 
methods where for each trial we produce simulated light curves 
taking the measured light curves as input and assuming Gaus- 
sian statistics. In Table ^j, we report the mean value of the fit- 
ted peaks and the 99.9% C.L.s for each pair of light curves. 
We have also estimated the cross-correlation function using the 
Lomb-Scargle method for computing the FFT of unevenly sam- 
pled data (Scargle 1989) and find essentially the same results. 

4. POWER SPECTRA 

We can combine our pointed X-ray observations, which span 
time scales from approximately 5 x 10 3 s to 3 x 10 5 s, with ob- 
servations from the All Sky Monitor (ASM) on board the RXTE 
spacecraft, which span time scales ranging from days to 3 years. 
We have used these data to construct an X-ray power spectrum 
that covers four orders of magnitude in Fourier frequency from 
10~ 8 -10~ 4 Hz. The result is presented in Fig. |[ A power spec- 
trum that covers a similarly broad range of time scales, but with 
much better statistics, has been obtained from pointed RXTE 
observations of the Seyfert 1 galaxy NGC 3516 (E delson & 
Nandra 1999). 

As the orbital time scales of the RXTE, ASCA, and EUVE 
spacecraft are approximately 5 ks, significant variability power 
(due to SAA passages and source occultation) appears on these 
time scales and at higher harmonics. In addition, for the EUVE 
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FIG. 3. — Power Spectral Densities constructed from RXTE-ASM, 
RXTE-PCA, and EUVE data (see text). Dashed lines show flat, f 1 , and 
f~ 2 overlaid on the power spectra. The dotted line show our estimate of 
the mean noise level of the ASM data. All PSDs are normalized such that 
integrating over positive Fourier frequencies yields the root mean square 
variability, normalized to the mean, for each particular light curve. 



spacecraft, we must integrate for approximately 5 ks in order 
to obtain good signal to noise. We therefore only consider 
Fourier frequencies < 10 -4 Hz. The PSDs are constructed from 
light curves with 5.6 ks time bins for the EUVE data, and 128 s 
time bins for the RXTE-PCA data. These light curves, how- 
ever, still contain a number of data gaps, therefore we use tech- 
niques for calculating power spectra o f unevenly sampled light 
curves (see Lomb 1976; Scargle 1982). The PSDs are then nor- 
malized so that integrating over positive frequencies yields the 
root mean square (rms) variability relative to the mean (see B el- 
loni & Hasinger 1990a; |Miyamoto et al. 1992| ). The resulting 
PSDs are then logarithmically binned over frequencies ranging 
from / — ► 1 .2/. For the EUVE PSD we calculate that the noise 
contribution to the rms variability is 1 1 .0 ± 0.6% of the mean, 
therefore we subtract this (assumed white) noise level from the 
EUVE PSD. The noise level of the RXTE PSD is negligible. 

Several results are immediately apparent in Fig. |[ First, we 
see that the EUVE light curve is more highly variable than the 
RXTE-PCA light curve, with an rms variability of 18.1 ± 1.4% 
compared to 7.4 ±0.6% over the time scales sampled by the 
PSDs. Both the RXTE-PCA and the EUVE PSDs are approxi- 
mately oc f~ 2 . The f~ 2 PSD seen for the RXTE-PCA light curve 
is comparable, in both shape and rms amplitude, to the high- 
frequency (> 3 Hz) PSDs seen for galactic black hole candi- 
dates ( GBHCs), such as Cygnus X-l, observed in t heir low/hard 
states ( |Miyamotoet al. 1992j ; |Nowak et al. 1999a| ). 

The low-frequency X-ray PSD is calculated via ASM data. 
The ASM provides light curves in three energy bands, 1.3- 



3.0keV, 3.0-5.0keV, and 5.0-12.2keV, ty pically consisting of 
several 90s measurements per day (see |b 



^evine et al. 1996 



Remillard & Levine 199% ^ochner & Remillard 1997| ). We 
have combined all three energy channels into a single channel, 
and then rebinned the data on 5 day time scales. Here we have 
filled data gaps (39 out of 233 points) with a linear interpola- 
tion, and the power spectrum is then constructed using standard 
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FFT techniques. We did not subtract the expected noise level 
calculated from the ASM error bars, as the derived PSD does 
not flatten at this level. This indicates that either the ASM er- 
ror bars are overestimated or that they do not follow gaussian 
statistics. 

The derived PSD is approximately fiat at low frequencies, 
shows a possible break at ss 7 x 10~ 7 Hz (w 200 day period), 
and is consistent with an /"' dependence at higher frequencies. 
The extent to which the PSD remains fiat at low frequencies is 
uncertain due to the limited time-span of the ASM observations 
compared to the putative 200 day break period. We note that 
Czerny, Schwarzenberg-Czerny, & Loska (1999) find a break in 
the optical power spectrum at a period <J 1000 days. The slope 
of the ASM-PSD at high frequencies is also uncertain due to 
the large noise level (shown as a dotted line in Fig. |]). The rms 
amplitude of the ASM lightcurve, including noise, is 30%. 
A Lomb-Scargle periodogram yields comparable results, and 
also shows a f ormally significant (via the methods of H ome & 
Bailunas 1986) peak at a 200 day period. This peak is coinci- 
dent with the break in the PSD, and furthermore the methods of 
Home & Bailunas are only strictly valid for a single sinusoidal 
period buried within pure counting noise. 

The overall X-ray PSD is again remarkably similar in both 
shape and rms amplitude to those seen in the low/hard states of 
GBH Cs. Low state GBHCs such as GX 339-4 (N owak et al. 
1999 b) or Cyg X-l flBelloni & Hasinger 1990a|; M iyamoto et al. 
1992; Nowak et al. 1999a| ) typically have an rms amplitude of 
30-40%, a flat PSD at low frequency with a break into an f~ x 
spectrum at 0.03-0.3 Hz, and a further break into an f~ 2 spec- 
trum at 1-10 Hz. The X-ray PSD we have found for NGC 5548 
is comparable to this, except scaled to approximately a factor of 
10 6 lower in frequency (although we do not resolve the break 
into the f~ 2 spectru m). A similar result was fo und for the X-ray 
PSD of NGC 3516 ( |Edelson & Nandra 1999| ). 

The factor of 10 6 is comparable to the expected ratio of the 
black hole masses between these two kinds of systems. Further- 
more, the 10 4 s time delays between the ASCA and hard RXTE 
bands discussed above are scaled by a similar factor compared 
to the 10~ 3 -10~ 2 s X-ray time lags seen at the 1-10 Hz break fre- 
quency in the low/hard state of GBHCs (see M iyamoto & Kita- 
moto 1989; |Nowak et al. 1999a| ). In the high-frequency regime 
in GBHCs, the varia bility coherence between X -ray bands be- 
gins to break down ( Vaughan & Nowak 1997|; N owak et al. 
1999a; |Nowak et al. 1999b[ ). In these regimes, we therefore do 
not expect a perfect linear correlation between energy bands. 
Carrying the analogy to GBHCs further, low/hard state GB- 
HCs exhibi t X-ray time lags of 0.03-0. 1 s at the lower break 
frequency (Miyamoto & Kitamoto 1989 ; Nowak et al. 1999a; 



Nowak et al. 1999t), hence one might expect OilOP s) time 



delays in the observed O(100 day) time scale variability. This 
is likely below the detection threshold of the RXTE- ASM, and 
we find no evidence for such delays comparing the lowest and 
highest ASM channels. Such time delays, however, if they were 
to exist, will likely be detectable by future all sky monitors with 
larger effective area and somewhat longer sampling time scales. 

5. X-RAY SPECTRAL FITTING 

In this section we use standard spectral fitting techniques to 
characterize the X-ray energy spectrum of NGC 5548. From 
previous X-ray studies, we expect the spectrum to be dominated 



by non-thermal power-law continuum emission. In addition, a 
relativistically broad iron emission line (from the surface of the 
accretion disk) and an associated reflection continuum will be 
superposed on this power-law continuum. These expectations 
guide our choice of spectral models (described below), and are 
confirmed in the sense that our spectral models describe the data 
very well. As we noted previously, we restrict our attention to 
data above 2 keV because of RDD effects; this removes the ne- 
cessity to model the warm absorber in this object. All spectral 
fitting has been performed using XSPEC V 10.00. 

In order to search for spectral variability, we examined X-ray 
spectra from ten distinct time intervals. For seven of these in- 
tervals, we have simultaneous ASCAIRXTE data: observations 
1, 3 and 4, as well as four sub-intervals of observation 2 (de- 
fined in Fig. [l]). The three additional intervals for which only 
RXTE observed occurred just prior to and after observation 3 
and during the period 16-18 August 1998 (cf. Tables | and||). 
This set of spectra allows us to probe both long (weekly) and 
short (daily) time scale spectral variability. The overlapping 
ASCA and RXTE data were pruned in order to ensure strict si- 
multaneity of the spectra from the two instruments. This was 
necessary since the observations for both satellites were inter- 
rupted either due to scheduling conflicts (RXTE) or telemetry 
failures (ASCA). However, since we are not examining spectra 
from data which are integrated over intervals of less than sev- 
eral tens of ks, we have neglected non-simultaneity on shorter 
time scales such as the orbital period. 

Since the count rates obtained from the RXTE-PCA observa- 
tions are much greater than that of the ASCA-SIS, we consider 
the RXTE-PCA data by itself initially. We then use those fits 
to guide our joint SIS/PCA fits, relying upon the greater spec- 
tral resolution of the ASCA detectors to constrain the iron line 
component better. 

5.1. The RXTE data 

We fit several models to these data. Our initial model is 
a simple power-law modified by Galactic neutral absorption 
(Nu = 1.7 x 10 20 cm" 2 ; z = 0). 1 This proves to be an unaccept- 
able fit in all ten cases, with residuals clearly indicating the need 
for an iron emission line and spectral hardening above lOkeV. 
Next, we add a narrow emission line to the power-law spectrum 
with a (rest-frame) energy fixed to that of the Ka line of cold 
iron (6.4 keV). Using the F-test, this leads to a significant im- 
provement (at greater than the 99% level) in the goodness-of-fit 
in all of the datasets. Adding more complexity to the fit, the 
iron line is allowed to have a broad Gaussian profile and arbi- 
trary energy. Again using the F-test, this leads to a significant 
improvement (at the 90% level) in the goodness-of-fit in two 
cases only (observations 1 and 2.2). 

Proceeding to more physical models, we replace the Gaus- 
sian emission line profile with the iron line profile expected 
from a relativistic accretion disk around a Schwarzschild black 
hole (Fabian et al. 1989). The inner edge of the line emit- 
ting region, the inclination of the disk, and the normalization 
of the line are free parameters in this fit. Hidden fixed param- 
eters are the outer edge of the line emitting region, r out , the 
line emissivity index (3 defined such that the local line emis- 
sivity is proportional to r~@, and the rest frame line energy 
which is fixed at 6.4 keV. The emissivity index (3 is fixed at 



'For all the fits to the RXTE-PCA data, we include a fixed power-law to represent the contamination by IE 1415.6+2557. The spectral indices used for this added 
power-law are shown in TableH. The flux normalizations at 1 keV shown in that table are reduced by a factor determined by the collimator response to IE 1415.6+2557 
for each pointing directed at NGC 5548. This factor is typically ~ 0.33. 
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= 2.5 which approximates the expected emissivity in a variety 
of physical scenarios. The outer radius r out is fixed at 1000r g , 
where r g = GM/c 2 is the gravitational radius of the central black 
hole. For our chosen value of (3, the fits are insensitive to the 
value of r out , provided that r out > 100r g . Although this is a more 
physical model, the RXTE-PCA data cannot distinguish it from 
the Gaussian line profile, i.e., there is no improvement in the 
goodness-of-fit in any of the datasets when the diskline profile 
is used. We do note, however, that both the broad Gaussian and 
the diskline fits are consistent with a substantial fraction of the 
iron emission originating from within ~ 20 r g . 

Our final models include the Compton reflection contin- 
uum which is implemented in XSPEC as the PEXRAV model 
(Magdziarz & Zdziarski 1995). We fix the e-folding energy of 
the underlying hard continuum at 120 keV as determined from 
fits by Magdziarz et al. (1998) to data from OSSE observations 
of NGC 5548, and we set the inclination of the planar cold re- 
flector at 30°. We also fix the chemical abundances to be So- 
lar. The relative normalization of this component, TZ, is left 
as a free parameter. For an isotropic primary X-ray source, 
TZ = Q./2-k, where f2 is the solid angle subtended by the cold 
reflecting medium as seen by an observer situated at the X-ray 
source. This leads to a large improvement in the goodness-of- 
fit for all of the observations, showing that the Compton con- 
tinuum (which produces the high-energy spectral hardening) is 
clearly detected. 

5.2. Joint ASCA/RXTE spectral fitting 

Four of our ASCA observations were scheduled to be con- 
temporaneous with RXTE observations. Joint spectral fitting 
of simultaneous data from ASCA-SIS and RXTE-PCA can be a 
powerful technique — it allows us to utilize both the spectral 
resolution of ASCA and the broad-band sensitivity of RXTE. 
The results of the joint ASCA-SIS/RXTE-PCA fits for our seven 
intervals are listed in Table ^. To facilitate direct comparison, 
we show the fits to the RXTE-PCA data alone which include a 
broad Gaussian emission line. The addition of the ASCA data 
does not affect the model fits greatly, although in some cases it 
allows the iron line width to be resolved. 

5.3. Observed correlations between spectral parameters 

Here we use the results of our spectral fitting of the ASCA 
and RXTE data to examine spectral variability. In order to en- 
sure that we have a set spectral parameters which have been 
uniformly analyzed, we consider only those parameters which 
have been obtained from fits to the RXTE data alone. Spectral 
parameters from the joint ASCA/RXTE fits yield similar results. 
Figure |] shows the photon index T, the iron line flux Ni, the 
iron line equivalent width W^a, and relative reflection normal- 
ization TZ as a function of the 2-10keV flux, F2-10, of the hard 
continuum component. For these plots and all of the fits to the 
spectral trends discussed in this section, we present uncertain- 
ties as 1-er error bars. 

5.3.1. The photon index 

There is a clear positive correlation between the photon index 
and the 2-10keV flux (Fig. |]a). Fitting a simple linear form to 
this correlation gives 

r = (1.63±0.04) + (0.030±0.006) [ — n Fl ~ W . * ) . (1) 

\10~ u ergs -1 cm -2 / 

This simple model does not strictly provide a formally ade- 
quate description of the spectral variation with flux, giving a 



goodness-of-fit of x 2 /dof = 16.1/8. However, we note that 
most of the contribution to this large value of \ 2 is due to the 
Aug 16-18 data point. Omitting this point, we find x 2 /dof = 
5.7/7 with the fit parameters being substantially unchanged. 
Nonetheless, the fact that the X-ray continuum gets softer when 
brighter is consistent with the trend observed in previous mon- 
itoring campaigns of this source (Magdziarz et al. 1998). How- 
ever, here we appear to find a somewhat weaker dependence of 
spectral index on flux. 

5.3.2. Iron line and reflection features 

Correlations of the iron line strength with the continuum flux 
are of great importance. In the simple X-ray reflection scenario, 
we would expect Wku to be constant provided that the light 
crossing time of the fluorescing region is much smaller than the 
time scale on which spectral variability is being probed. If the 
light crossing time of the emission line region is much greater 
than the time scale being probed, a constant iron line flux, Ni, 
would be expected. When these two time scales are compara- 
ble, reverberation effects come into play and both Ni and Wx a 
will be seen to vary. 

Like Wku, the parameter 1Z should measure the amount of 
X-ray reflection relative to the direct continuum. Indeed, Wx a 
should be proportional to 1Z provided the following conditions 
are satisfied: 

1. The Compton reflection continuum is not starting to 
dominate the observed continuum at the iron line ener- 
gies. In practice, this condition implies 1Z < 2-3. 

2. The ionization state of the illuminated regions of the re- 
flecting medium is fixed. 

3. The primary continuum has a fixed energy spectrum. Of 
course, we have just shown that this is not strictly the case 
and that the photon index of the primary continuum emis- 
sion changes by AT ~ 0.2 during our campaign in a man- 
ner that is well correlated with the 2-10 keV flux. How- 
ever, the Monte Carlo simulations of George & Fabian 
(1991) show that such small photon index changes only 
have a small (less than 10%) effect on the iron line equiv- 
alent width. 

Apart from inclination dependences (which are usually weak) 
and ionization gradients, the geometry of either the reflecting 
medium or the X-ray source does not affect the expected pro- 
portionality of 1Z and Wfc a . 

Our spectral results hint at a more complex picture. Inspec- 
tion of Fig. Qb shows that the iron line flux is essentially con- 
stant. This is borne out by the linear (dotted line) and con- 
stant (solid) models which we have fit and which describe these 
data very well. Apart from effects introduced by the subtle 
changes in photon index, this implies that Wk u is inversely pro- 
portional to the flux (Fig. ^c). More formally, fitting the 10 
data points in the Wku-F2-\o plot with a constant model results 
in x 2 /dof = 20.8/9 (dotted curve), unacceptable at the 97% 
level. Assuming, instead, a power-law relationship between 
these two parameters such that Wko °c ^2-10 results in a dra- 
matic improvement in the goodness-of-fit with x 2 /dof =7.9/8 
and a = -0.9 ±0.4 (solid curve in Fig. Qc). Thus, these results 
are consistent with a constant line flux and an equivalent width 
which is inversely proportional to the continuum flux. 
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FIG. 4. — (a) Photon index T, (b) iron line flux Ni, (c) iron line equivalent width Wx a , and (d) relative reflection normalization 1Z as functions of the 2-10 keV hard 
continuum flux for our ten time intervals. The iron line flux is in units of 1CT 5 ergs cm -2 s , and the 2-10 keV hard continuum flux is in units of 10 _1 1 ergs s cm -2 . 
The dotted curves are the expected relationships in simple reflection models, while the solid curves are the best fits assuming various functional forms (see text). 
One-sigma error bars are shown. 



Since the relative normalization of the Compton reflection 
continuum 1Z should be proportional to Wk u under a fairly gen- 
eral set of conditions, we would expect that 1Z oc F^q- Fitting 
this relationship to the data depicted in Fig. |]c gives x 2 /dof = 
14.5/9 (dotted curve). This is formally inconsistent with the 
data at the 91% level. Using a more general power-law form 
1Z oc f 2 -io leads to a vastly improved fit of \ 2 / dof = 3 .4/8 with 
a = 0. 14 ± 0.60 (solid curve). This is consistent with a constant 
relative reflection normalization and implies that Wxa and 7Z do 
not exhibit the proportionality which is normally expected. In 
order to examine this explicitly, we plot Wku vs 1Z in Fig. |^. The 
dotted line is the linear relationship expected for a cold reflec- 
tor with solar abundances for which Wxa = 150eV when 1Z = 1 
(George & Fabian 1991). This model yields a x 2 /dof = 73/10 
and is definitely in conflict with the data. However, if we fit for 
the proportionality constant, we obtain Wxa = 259 eV at 1Z = 1 
and x 2 /dof = 14.4/9 which is still formally inconsistent with 
the data at the 89% level. 



6. IMPLICATIONS FOR DISK/CORONA MODELS 

If we assume that the EUVE light curve is representative of 
the seed photons for the higher energy X-rays, we can make a 
simple estimate of the size of the scattering region based upon 
the delays seen. With each scattering, a seed photon picks up 
a fractional energy AE/E ~ 4kj,T /m e c 2 , where T is the tem- 
perature of the electrons in the corona. The energy of a Comp- 
tonized photon will be 

/ 4k b T\" 
E~E 1 + , (2) 

\ m e c L ) 

where Eq is the initial seed photon energy and n is the num- 
ber of scatterings (Rybicki & Lightman 1979). The time delay 
of this photon with respect to the seed photon source will be 
roughly proportional to the number of scatterings, t ~ ntQ. Here 
to ~ h/c where lj is the mean free path for Thomson scatter- 
ing (or the size of the corona if it is optically thin). We take 
the effective photon energy in the 0.5-1 keV ASCA data to be 
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FIG. 5. — Iron line equivalent width WKa versus the relative Compton 
reflection normalization 1Z. The dotted line is the expected linear relation- 
ship assuming solar abundances for the cold reflector (George & Fabian 
1991). The solid line is the best-fit proportionality relationship which is 
still in conflict with the data at the 89% level. 

E a = 0.78 keV and the effective photon energy of the 2-20 keV 
RXTE-PCA data to be E pc . d = 5.4keV. Both of these energies 
are determined using the instrument sensitivities in the rele- 
vant bands and assuming a photon spectral index of 1.9 and 
a neutral absorbing column of 1.7 x 10 20 cm 2 . Using a value of 
kt,T ~ 50keV which is found from fits to the OSSE spectrum of 
NGC 5548 (Magdziarz et al. 1998) and the inferred lags of the 
RXTE-PCA and 0.5-1 keV ASCA-SIS light curves relative to 
that of the EUVE, we solve for the characteristic time between 
scatterings and find to ~ 3.8ksec. 2 For a homogeneous corona 
with Thomson depth of r ~ 2 this implies a coronal size scale 
of ~ 2 x 10 14 cm, which is of order 10r g for a 10 8 M Q black 
hole. 

This interpretation for the relationship between the EUV and 
higher energy X-ray flux is not without its complications. First 
of all, as we note above, the fractional size of the step in the 
EUV is larger than the steps in the higher energy bands. Sec- 
ond, the amplitude of the EUVE PSD differs from that of the 
RXTE-PCA (Fig. ||). Therefore a strictly linear relationship be- 
tween the fluxes in these bands is not possible. One possibility 
is that the EUVE observations are sampling the Wien tail of a 
black-body spectrum so that the observed variations are more 
pronounced and are not necessarily characteristic of the bulk 
of the seed photons. A more provocative explanation would be 
that the increase in soft photon flux produces excess cooling in 
the corona, changing some combination of the corona temper- 
ature, size, or optical depth. The softening of the underlying 
continuum with increased flux (Fig. ^) does in fact show that 
properties of the corona must be changing. 

The iron line fluxes and reflection fractions also have signif- 
icant implications for the corona/disk geometry. As we noted 
above, both the reflection fraction and the iron line equivalent 
width are substantially smaller than has been previously seen 

2r 



for this object as well as for other type 1 Seyferts (e.g., MCG-6- 
30-15, Lee et al. 1998; and IC 4329A, Madejski et al. 1995), 
for which values of 1Z ~ 1 and Wxa > 160 eV are more typical. 
The lower values we find for NGC 5548 would be consistent 
with a substantially smaller covering factor for the cold mate- 
rial. This could occur, for example, if the Comptonizing region 
lay inside an inner disk radius as has been posited for models 
of Cygnus X-l (Dove et al. 1997; Gierliriski et al. 1997). Alter- 
natively, the disk could extend beneath the corona but may be 
completely ionized preventing this material from contributing 
to either the line fluorescence or the reflection hump. Mildly 
relativistic motion of the Comptonizing medium away from 
the disk could also account for the lower reflection fractions 
and iron line equivalent widths (Reynolds & Fabian 1997; Be- 
loborodov 1999). 

However, a lack of proportionality between the iron line 
equivalent width and reflection fraction (Fig. ||) cannot be ac- 
counted for simply by the relative geometry of the disk and 
corona. A constant iron line flux would suggest that a sub- 
stantial portion of that emission is produced very far from the 
Comptonizing X-ray source, perhaps in the outer disk or ob- 
scuring torus, so that the line emission does not vary on time 
scales as short as that of the underlying continuum. We would 
then expect that the reflection hump is likewise produced in this 
same distant material so that its absolute normalization should 
also be constant. Unfortunately, there is the additional com- 
plication that the Gaussian and diskline fits to the fluorescent 
iron line indicate that the line is broad and redshifted. The line 
shapes are consistent with a significant fraction of this emission 
originating from an accretion disk within ~ 20 r g . Furthermore, 
spectral fits to the ASCA-SIS data which include a narrow line 
at 6.4 keV along with the diskline model yield a narrow com- 
ponent that can contribute only about 15% to the total line flux. 
Therefore, based on the spectral shapes of the iron line, a sig- 
nificant amount of reprocessing by distant cold material is un- 
likely. 

Wherever the cold reprocessor is located, these results appear 
to be somewhat in conflict with models in which the spectral 
variability of the underlying continuum is linked to the relative 
reflection normalization 1Z (Magdziarz et al. 1998; Zdziarski 
et al. 1999). In Fig. |[ we plot the relative Compton reflection 
normalization versus the photon spectral index of the under- 
lying power-law (filled squares, solid contours). As we note 
above, these data are formally consistent with a constant value 
of 1Z, but there does seem to be a slight linear trend in the same 
sense predicted by the models. For direct comparison, we also 
plot in Fig. ^ the spectral parameters found by Zdziarski et al. 
(1999; see also Magdziarz et al. 1998) which are derived from 
Ginga observations. In this figure, we show the l-cr error con- 
tours for both sets of data. The spectral parameters we find 
appear to occupy a different region of the 1Z-T plane than that 
found by these previous observations. However, consideration 
of systematic uncertainties may alter this conclusion. Ginga- 
LAC observations of the Crab have yielded a photon index of 
T = 2.08 ±0.03 (Turner et al. 1989) which agrees with early 
rocket measurements by Toor & Seward (1974). In contrast, 
power-law fits to RXTE-PCA observations of the Crab yield a 
significantly softer index of T = 2.187 (Wilms et al. 1999). If 
this difference of AT ~ 0.1 is due to a systematic shift in mea- 
sured power-law spectral index, then one may be tempted to 
shift each of our values of T downward by this amount, making 



"One can also use the direct lag measured between the A5CA-SIS and RXTE-PCA light curves of 5 ks. From that lag, one finds fo ~ 0.9 ks. 
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FIG. 6. — The relative Compton reflection normalization 1Z versus pho- 
ton spectral index T for our fits to the RXTE-VCA data (filled squares, 
solid error contours) and for fits to Ginga data as presented in Zdziarski 
et al. (1999; unfilled squares, dotted contours). Although the RXTE-PCA 
parameters are formally consistent with a constant value of 1Z, a linear 
trend may be present. However, the dependence on T found for these data 
is weaker that that found by Zdziarski et al. (1999), and these parameters 
appear to occupy a different part of the 1Z-V plane. One-sigma contours 
are shown. 

the two sets of parameters consistent. Unfortunately, the ori- 
entation of the error contours clearly indicates a correlation be- 
tween fitted values of 1Z and T. This seems to be an intrinsic 
feature of the model. If so, any systematic shift of T probably 
also implies a bias in the fitted value of 1Z as well. It is there- 
fore not entirely clear whether these two datasets can be made 
consistent even if one considers the possible systematic errors. 

As a final caveat, we note that the strength of the reflec- 
tion component can be masked by features in the hard con- 
tinuum which need not be well-represented by a power-law. 
Stern et al. (1995) point out that for Comptonizing regions 
which are localized and which are illuminated anisotropically 
the resulting hard continuum will have spectral breaks just 
above the peak energies of the twice scattered photons (see also 
Haardt 1993). For sufficiently high values of the compactness 
(= LdissOY / hm e c 3 where Ldiss is the power heating the corona 
and h is the coronal size scale), these breaks may occur at ener- 
gies such that the reflection component will be underestimated. 
In addition, as the compactness increases, the spectrum will 
soften and the break will move to lower energies. Hence, any 
variability in the absolute normalization of the reflection com- 
ponent which is inferred from PEXRAV-type spectral fits could 
be significantly affected by changes in the compactness. This 
sort of uncertainty can only be resolved by more sophisticated 
spectral models. 



7. SUMMARY AND CONCLUSIONS 

Here we summarize the main results of this work: 

1. The EUV leads the harder X-rays, rejecting scenarios in 
which the soft X-ray component is produced by repro- 
cessing of harder X-rays. In this respect, these hard X-ray 
lags are reminiscent of the X-ray lags seen in NGC 7469 
(Nandraetal. 1998). 

2. However, in contrast to NGC 7469, the shorter lags we 
find for NGC 5548 are consistent with Compton dif- 
fusion time scales for a relatively small corona of size 
~10 14 cm. 

3. The power spectrum of the X-ray variability shows a 
break at ~ 200 day time scales. The PSD is similar to 
that seen in NGC 3516; and if one assumes a 10 8 M Q 
black hole for NGC 5548, this scales appropriately with 
mass compared to PSDs seen for the low/hard states of 
GBHCs such as Cyg X-l or GX 339-4. 

4. The iron line equivalent width and relative reflection nor- 
malizations are smaller than those which are typically ob- 
served for type 1 Seyferts. 

5. The broad, redshifted iron line profiles, the variability of 
the iron line equivalent width, and the apparent constancy 
of the relative reflection normalization are difficult to rec- 
oncile in the context of simple reflection models. 

The last point, while not definitively constrained by the data, 
may pose a serious challenge for models of reflection by cold 
material and will only be resolvable with further detailed tem- 
poral and spectral observations. More sensitive studies with 
future instruments such as XMM would be extremely useful to 
verify this result. 
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EUVE/ASCA/RXTE Observations of NGC 5548 



Obs. 


Date 


Exp. time (ks) 


mean count rate (10 2 cps) 


Spectral Density at 76A(^Jy) 




02 Jun 


8.8 


18.0±0.6 


345 




09 Jun 


11.5 


8.4 ±0.4 


161 


2 


18 Jun 


139.8 


9.3±0.1 


178 


3 


01 Jul 


21.3 


2.1 ±0.2 


40 



Table 1 
euve-ds observations 



Obs. 


Date 


good SIS exp. 


good GIS exp. 


SIS0 count 


GIS2 count 


2- 


-lOkeV flux 






time (ks) 


time (ks) 


rate (cps) 


rate (cps) 


(io- 


11 erg cm" 2 s" 1 ) 


1 


15 Jun 1998 


21.2 


25.3 


2.61 


1.56 




5.9 


2 


20 Jun 1998 


111.1 


106.3 


3.07 


1.71 




6.9 


3 


01 Jul 1998 


9.2 


10.6 


1.66 


0.97 




4.2 


4 


07 Jul 1998 


14.7 


16.9 


1.73 


1.08 




4.3 




19 Jan 1999 


8.2 


9.0 


1.84 


1.09 




4.7 



Table 2 

Basic parameters of our ASCA observations of NGC 5548. The quoted fluxes are those observed (i.e. subject to the 
total line-of-sight absorption). A simple power-law model was used to measure these fluxes. 



Obs. 


Date(s) 


target 


Good Time (ks) 


PCA/PCUs 0-2 


HEXTE rate 










rate (cps) 


(cps) 


1 


Jun 15-16 


NGC 5548 


21.24 


25.9 


1.0 






IE 1415.6+2557 


8.82 


4.8 




2 


Jun 19-24 


NGC 5548 


110.28 


28.7 


0.7 






IE 1415.6+2557 


7.56 


4.4 








NGC 5548 scan 


3.48 






3 


Jun 29-Jul 2 


NGC 5548 


41.70 


19.3 


0.6 






IE 1415.6+2557 


7.44 


5.2 




4 


Jul 7 


NGC 5548 


24.60 


19.3 


1.0 






IE 1415.6+2557 


3.60 


4.2 






Aug 16-18 


NGC 5548 


44.92 


23.8 


0.3 






IE 1415.6+2557 


7.74 


7.3 








Table 3 










RXTE observations 










Instruments 


fi ag (ks) 99.9% C.L. (ks) 








EUVEvsASCA 


13.1 -3.9 


29.5 








EUVE vs RXTE PCA 


35.3 23.1 


53.1 








ASCA vs RXTE PCA 


5.0 1.9 


8.0 





Table 4 

Lags obtained using the ZDCF for the three pair-wise combinations of instruments during the 18-23 June 1998 

observation period. 



Obs. 


PCA/PCUs 0-2 rate (cps) 


r 


A (10" J ) 


1 


4.8±0.2 


2.3±0.1 


9.8(-1.3/ + 3.8) 


2a 


3.9±0.3 


2.4 ±0.2 


11.6(-2.8/ + 3.8) 


2b 


4.8±0.3 


2.2 ±0.2 


8.8(-2.1/ + 2.8) 


3a 


5.4±0.2 


2.3±0.1 


ll.l(-2.0/ + 2.4) 


3b 


4.9±0.2 


2.2±0.1 


9.3(-1.6/ + 2.0) 


4 


4.2±0.3 


2.3±0.1 


11.6(-2.3/ + 2.9) 


Aug 16-18 


7.3±0.2 


2.1 ±0.1 


10.8(-1.0/+ 1.0) 



Table 5 

Properties of IE 1415.6+2557. 
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Obs. 


-F2-10 


r 


n 


Gaussian 




X 7 /dof 






SIS PCA 




E (keV) a (keV) 


N, 





0.08 
+0.20 



1 

2.1 
2.2 
2.3 
2.4 

Jun 29 

3 

Jul 2 

4 

Aug 16-18 



8.0 
7.9 
7.8 
7.7 
9.8 
9.7 
9.5 
9.4 
9.1 
8.8 
6.9 
5.6 
5.3 
5.4 
5.8 
5.7 
7.2 



1.86 ±0.03 0.34 



1.89 ±0.03 



1 op+0.05 
1 Q4+0.03 



1.93 ±0.02 



1 Qc+0.06 
l - yJ -0.04 



1 70+O.O6 



1 77+O.O5 
i - / '-0.01 



1.86 ±0.03 

1.86±0.03 

1.93 ±0.02 

1.91 ±0.02 

1 90 +002 

1 90 +002 
I 91+0.04 

L9l|l 
1-80!^ 

1 09+0.04 

i ' oz '-0.02 

1.80 ±0.04 

1.77 ±0.03 
1 7C+0.03 
i - ' J -o.oi 
1 00+0.01 

1 -°°-0.02 



0-4o: . 07 

0.385" 

0-40tl* 

0.46±0.12 

41+0-16 
u -^ 1 -0.10 

0.3733 

0.40!^ 
46 +018 
49+0-30 

29 +0 29 

n to+0.18 

62 +0 33 

n o 1 +0.40 
u - 31 -0.30 

33+ 018 

U - JJ -0.14 

32+ 20 
43+ 05 



6 1 5+ ulu 
6 1S+ 013 
6 15+ 011 
6 17+010 

6.10±0.12 
6.11±0.12 

(T 1O+0.13 

6 17+°- 15 

' —0.06 

6 32 +0 - 22 
6 25 +0 - 25 

°- ZJ -0.20 



<0.41 
3S+ 017 
0.33^ : [f 

0.30^8:11 

0.44!0;{3 
3S+ - 26 

u ' JJ -0.09 
U-^-O.ll 

41+0.12 
u -^ 1 -o.io 
<0.74 
46 +0 ' 30 

u - w -0.23 



6.21 ±0.12 <0.36 



6 20+ 009 
6.16±0.14 
6 21 +018 

°' Z1 -0.14 

6 17 +009 

D1 '-Q.1Q 

6.08^ 



0.38^8:11 

0.44!°;|I 

< 1.11 

39 +a 12 
U - Jy -0.10 

29 +014 

u ' zy -0.12 

<0.17 



8-i^: 

10.6+11 
8.1!!;? 

Q 1+2.0 

9.6!l| 

8.1 ± 1.3 

7 5 +5 - 7 
-3-1 

6 2 +34 

8.8! 2 | 
q 1+2.0 
1.8 

7.3 ±1.4 

8-2l| 
9.8+ 20 , 

7 +o:9 
'• u -i.o 

7 5+ L1 



47.4/40 
382/381 
38.4/40 
327/324 
36.0/40 
501/409 
29.8/40 
423/416 
32.2/40 
241/286 
26.0/40 
28.0/40 
202/278 
26.3/40 
47.0/40 
287/332 
29.1/40 



Obs. 


F2-10 




r 


K 


DISKLINE 




xVdof 






SIS 


PCA 




n n i (deg) 


Ni 





1 

2.1 

2.2 

2.3 

2.4 

3 

4 



7.9 
7.7 
9.6 
9.4 
8.8 
5.4 
5.7 



1.89 ±0.03 
1.88 ±0.04 

1.95 ±0.02 
1.93 ±0.02 

1.96 ±0.05 

1 70+O.O6 
1 - /8 -0.04 



1.86 ±0.03 
1.85 ±0.03 
1.91 ±0.02 

1.90 ±0.02 

1.91 ±0.04 
1 o 9 +o.o3 

l- 8Z -0.04 



-0.18 



1.77 ±0.04 1.75 ±0.03 



0.39 

39+ o:l9 

0-42!";!° 
0.38i} 2 
5O +0 - 28 
<S9+°- 27 
32 +018 

U -- 3Z -0.16 



3.5 



9.5^ 
10.3!f 2 

7 6 +3 - 7 

8 9+ 67 
12. 1 +77 



< 



6.1 
16.4 



1 1 3 +70 

1 L - 3 -3.3 



<32 
<28 
<29 

26l;f 

4138 
<37 

<25 



6.8! I 
6-8!ii 

6.1^ 

7 o+1.6 
'■ J -1.7 



383/381 
327/324 
503/409 
420/416 
242/286 
203/278 
288/332 



Table 6 

The parameters obtained from spectral fits to the RXTE and ASCA data. The spectral models consist of Galactic 

ABSORPTION, A FLUORESCENT IRON EMISSION LINE, AN UNDERLYING CUT-OFF POWER-LAW AND A COMPTON REFLECTION 
COMPONENT. FOR FITS IN THE UPPER PART OF THE TABLE, A BROAD GAUSSIAN FUNCTION IS USED TO MODEL THE IRON LINE; IN THE 
LOWER PART, A DISKLINE MODEL IS USED. FOR THE JOINT ASCA-SIS/RXTE-PCA FITS, WE REPORT SEPARATE SPECTRAL INDICES FOR 
THE UNDERLYING POWER-LAW FOR SIS AND PCA. THE 2-10 KEV CONTINUUM FLUX IS REPORTED FOR THE PCA FITS ONLY AND IS IN 
UNITS OF 10"" ERGSCM~ 2 S~'. THE IRON LINE FLUX, Ni, IS IN UNITS OF 10 5 ERGS CM 2 S 1 , AND THE INNER RADIUS OF THE DISKLINE 

MODEL, r m IS IN UNITS OF r g = GM/c 2 . 



